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Abstract —In this article we study synchronization of Ku¬ 
ramoto oscillators with heterogeneous frequencies, and where 
underlying topology is a graph of diameter two. When the 
coupling strengths between every two connected oscillators are the 
same, we find an analytic condition that guarantees an existence 
of a Positively Invariant Set (PIS) and demonstrate that existence 
of a PIS suffices for frequency synchronization. For graphs of 
diameter two, this synchronization condition is significantly better 
than existing general conditions for an arbitrary topology. If 
the coupling strengths can be different for different pairs of 
connected oscillators, we formulate an optimization problem that 
finds sufficient for synchronization coupling strengths such that 
their sum is minimal. 

I. Introduction 

Synchronization of coupled oscillators is an important topic 
of research for scientists from different areas including neuro¬ 
science [1], (28), (29), physics [3), (23| , mathematics (14), and 
engineering (9), |20| , Kuramoto model [151 of coupled 
oscillators, despite its seeming simplicity, demonstrates a quite 
rich dynamic behavior and has become a canonical model for 
studying synchronization. 

The two main features that describe the behavior of a 
system of coupled oscillators are the coupling function and 
the interconnection topology. In the case of the Kuramoto 
model a trigonometric sin() is used as the coupling function; a 
broader class of the coupling functions, however, has also been 
discussed d), CD. G3- The most popular assumption 
on the interconnection topology is that all oscillators are con¬ 
nected to each other, which corresponds to a fully connected 
graph or a graph of diameter one 0,0- A much more general 
approach is to study the systems of oscillators with an arbitrary 
underlying topology 0, §, (12), (T4), (24)- 

Several additional assumptions can be made to make 
analysis of the Kuramoto model more tractable. First, one 
may consider a limit case when the model contains infinite 
number of oscillators gd> gd 1 16), [26). Second, it can be 
assumed that all oscillators have equal intrinsic frequencies, 
and therefore form a gradient system of homogeneous oscilla¬ 
tors fl8) . Alternatively, as we do in this article, one may let the 
frequencies to take distinct values and thus analyze a system 
of heterogeneous oscillators ©, ©, 0,d), O, (27). 
Finally, the coupling strengths can be equal for all pairs of 
connected oscillators, or are allowed to take different values 
for different connections. 


In this paper, we consider a system of finite number of 
heterogeneous Kuramoto oscillators in which the underlying 
topology is a graph of diameter two, a natural step to further 
generalization of the complete graph (diameter one) case. First, 
we consider the case when the coupling strength is the same for 
all pairs of connected oscillators and formulate an analytic con¬ 
dition that guarantees boundedness of the trajectories, which in 
our case also implies synchronization. While there exist more 
general synchronization conditions 0, 0, 0, (T?), (24) that 
are applicable to the systems with an arbitrary topology, they 
are significantly more restrictive (when applied to the diameter 
two graphs) compared to our analytic condition. We provide 
simulation results that illustrate the improvement over existing 
results for the graphs of diameter two. Second, when the 
coupling strengths are allowed to be different for different pairs 
of interconnected oscillators, we formulate an optimization 
problem that finds the coupling strengths such that their sum 
is minimal while synchronization is preserved. 

The rest of the paper is organized as follows: in Section II 
we describe the problem setup as well as the main challenge 
for guaranteeing synchronization, i.e. showing boundedness 
of trajectories. This challenge is addressed in Section III- 
A, where we show a general, yet hard to check, condition 
for synchronization (Proposition 1). This condition is made 
tractable in Section III-B for the case of equal coupling 
strengths. Further, in Section III-C we present an optimization 
approach to study the case when the coupling strengths can 
be different for different pairs of connected oscillators. We 
illustrate our findings using simulations in Section IV and 
conclude in Section V. 

II. Problem Formulation 

In this article we study a system of Kuramoto oscillators in 
which each oscillator is described by the following equation: 

<Pi=uJi+ Y — • sin {(j)j - <j>i), (1) 

“ n 

jeNi 

where Ni is a set of oscillators connected to oscillator i, 
i.e. the set of its neighbors, Kij is the coupling strength 
between oscillators i and j, and n is the total number of 
oscillators in the system. The coupling strength is symmetric 
(K^ = Kji, Vi, j), and can be the same for all connections 
as assumed in Section III-B, or can be different for different 


pairs of connected oscillators as in Section III-C. We also 
assume that the intrinsic frequencies of oscillators uji are 
heterogeneous, which implies that they are not necessary equal. 
Frequencies, however, do not change their values with time, 
so each Wi is a constant. 

In this article we study frequency synchronization of the 
system 0. System q achieves synchronization if = 

<p 2 (t) = • • • = <j) n (t) as t -A oo. We will denote the common 
phase velocity by </>. 

This common phase velocity <j) is an average intrinsic 
frequency of the oscillators: 

n 

E 

4 >= . 

n 

Indeed, when = </> 2 = • • • = m <j> n , the sum of all the 
equations of ([TjMs: (</>i + 02 + • • • 4>n) = tui H - • • • M n -> 

because each sin — 4>i) is added to sin(<^ — <fij) 
and gives zero. 

We now define the deviations of the frequencies for all 

n 

i = 1,..., n by Oji = uji — Oj, where Oj = - is the 

n k =1 

average intrinsic frequency, and study the following system: 

(f>i=Ui+ Y] — • sinO j - (2) 

jeNi 

Each limit cycle of system ([I]) is an equilibrium of ([2]). 
Therefore, we will focus on finding conditions when system 

0 synchronizes, i.e. when 0* s= 0 Vi * 1,_, n. Due to the 

rotational invariance of system ([2]), and since the phase average 
remains the same (</>i H-b 4> n = 0), we can assume without 

n 

loss of generality that = 0 for all t > 0, where are 

i=l 

the trajectories of system ([2]). 

We will show frequency synchronization of system ([2]) by 
providing a Lyapunov function and using LaSalle’s Invariance 
Theorem G3- When the oscillators are homogeneous, all the 
intrinsic frequencies are equal, i.e. deviations uj% = • • • = 
O n = 0, and the following Lyapunov function can be used: 

V 0 {(j>) = - ^2 • cos (<t>i - fa), 

ijeE,i<j 

where 0 E M n and E is the edge set of a given graph. It can 
be verified that: 

n 

^o($ = -X>?<° 

Since function Vo(0) is well-defined on a n-dimensional torus 
T n which is compact, applying the LaSalle’s Invariance The¬ 
orem (on T n ) guarantees synchronization of the oscillators. 

When the intrinsic frequencies are not equal, we have a 
system of heterogeneous oscillators, and we still can provide 
a potential function for this case: 

v{$) = - - E ■ ( cos (^ - ^))- 

k= 1 ij£E,i<j 


We can check that the derivative of this function is also non¬ 
positive and is equal to zero only at equilibrium, i.e. when the 
frequencies are synchronized. 

The key problem here is that function is not bounded 
from below and cannot be defined on T n . Therefore, we are 
not able to apply directly the LaSalle’s Invariance Theorem. 
However, if we show that the trajectories <p E M n of 0 are 
bounded, then the function V(4>) is bounded as well, hence 
synchronization follows. 

One of the techniques for showing boundedness of the 
trajectories is to find a bounded Positively Invariant Set (PIS) 
for the oscillators’ phases. The goal of this article is to show 
that when some conditions are met, a PIS exists, and if the 
initial phases are in this PIS, then the trajectories will be 
bounded and, therefore, system 0 will achieve frequency 
synchronization. 

III. Main Results 

This section is organized as follows: we first introduce the 
notations used in this article and provide a general synchro¬ 
nization condition in Proposition 1. We also demonstrate by 
means of an example that existence of an equilibrium does not 
guarantee that system 0 achieves frequency synchronization 
for all initial phase values. In Subsection B we provide an 
analytic synchronization condition for system 0 with equal 
coupling strengths. In Subsection C we study a more general 
case when the coupling strengths can be different for different 
edges. 

A. Preliminary Results 

Let G = (V, E) be an undirected graph with vertex set V 
and edge set E that defines the topology of the system 0. 
Distance between vertices i and j is defined as a number of 
edges in the shortest path between i and j , where the length 
of a path is defines as the number of edges in it. Diameter of 
a graph is defined as the maximum distance between its two 
vertices. All the results presented in this article are formulated 
for the graphs of diameter two. 

We denote by A the symmetric adjacency matrix of a graph 
G, and define for each pair of vertices i, j constant P^ : 

Pij — a i * + 2 Aij , (4) 

where ai and aj are the i th and j th rows of matrix A and A^ 
is the ( ij) th element of matrix A. The dot product ai • a? 7 is 
equal to the number of common neighbors of vertices i and 
j, and Aij = 1 if and only if there is an edge between i 
and j in E. For example, if i and j are connected and have 
3 common neighbors, then P l3 = 5. Since diameter of the 
graphs considered in this article is less than three, > 1 for 
all pairs of vertices i, j. 

We denote the maximum and minimum phase values at 
time t by </4 aa , = max</>* and <f>\ 3in = min$, where <p\ is a 

^ l 

phase of oscillator i at time t. Let D t be defined as a maximum 
phase difference between two oscillators at time t (t > 0), i.e. 

Dt tf)max 4*mim 



then $min < <t>\ < 4>max (Vi = 1, ■ ■ ■ ,n). In other words, 
each phase lies between the minimal and maximal phases fa min 
and fa max . Maximum initial (at time t = 0) pairwise phase 
difference is denoted by Do: 

Do 0 max truin' 

If we can show that the maximum phase difference is always 
bounded, i.e. if D t < D V t > 0, where D is a constant 
satisfying D 0 < D < oo, then the trajectories will be also 
bounded since the phase average remains the same (for system 

([2]): fa 2 -|-0n = 0)- The PIS, therefore, is defined through 

the maximum phase difference that is bounded by the value 
of D: 

n 

PIS = {0 E M n : max \fa — fa\ < D, fa = 0}, (5) 

hj .. 

i=i 

which is obviously a compact. 

We now formulate a general sufficient condition that guar¬ 
antees that the maximum phase difference is always bounded 
by a constant D and thus the trajectories are also bounded. 

Proposition 1 If D is a constant satisfying Do < D < oo, 
and for all times t > 0 such that D t = fa max — 0^ in = D, 
the following condition is satisfied: 

4>k ~ 4>\ = - ui 

- *S2 ~ ■ sin (^fe -$)- ^2 —■ sin (^i - <t>\) < 

ieN k 71 jeN, 71 

(6) 

for every two oscillators k and l such that fa, = 0^ aa , and 
fa = famin, then the maximum phase difference is bounded by 
D, i.e. D t < D for all t > 0, trajectories of system 0 are 
bounded, and system 0 achieves frequency synchronization. 

Proof: Condition <01 says that when the maximum phase 
difference achieves value D, it can not grow anymore and 
thus does not exceed D. This implies that the trajectories of 
system 0 are bounded in W 1 since the phase average is always 
equal to zero. Further, function V (0) is well-defined in M n and 
we can apply LaSalle’s Invariance Theorem to guarantee that 
each solution of ([2]) approaches the nonempty set {V = 0} = 
{fa = 0.1 < i < n}, and system 0 achieves frequency 
synchronization. ■ 

It is possible that when 0^ aa , — fa m i n = D, several 
oscillators have phase values equal to 0^ aa . or fa min . In 
this case condition should be satisfied for each pair of 

oscillators with a phase difference equal to D. 

Condition 0 is very general by itself and cannot be 
directly applied to ensure boundedness of the trajectories and 
frequency synchronization of a given system. In the next two 
subsections we derive two conditions that can be easily verified 
for each given system and guarantee that condition ([6]) of 
Proposition 1 is satisfied. In particular, in Subsection B we 
derive an analytic condition for the case of equal coupling 
strengths, and in Subsection C we formulate an optimization 
problem for the case of non-equal coupling strengths. 

An alternative line of works ©- m focuses on results 
that guarantee existence of a locally stable equilibrium man¬ 
ifold for system 0. These local results, however, cannot 


1 
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Fig. 1: Topology in Example 1: a star graph with 3 nodes 




Fig. 2: Trajectories fat), t > 0 (left), and V(fa) (right) in 
Example 1 


guarantee synchronization for any given values of the initial 
phases (different from the equilibrium phases). We finish this 
subsection with an example that demonstrates that existence 
of a locally stable equilibrium for system ([2]) does not imply 
synchronization of this system for all possible initial phases. 
Therefore, existence of an equilibrium is not a sufficient 
condition of synchronization for all initial phases. 

Example 1 In this example three oscillators are connected 
as shown on Fig. |T] i.e. they form a star graph with three nodes. 
We assume that udi = 2 — e, = ^3 = — 1 + f, where e is 
a small positive constant, and all coupling strengths are equal: 
K 12 = K \3 = 3. It is easy to verify that this system possesses 
a locally stable equilibrium: fa — | sin -1 (l — §), fa = 03 = 
— | sin _1 (l — |). However, there are initial phases 0?, fa and 
03 for which the system does not achieve synchronization. On 
Fig. [ 2 ] the behavior of oscillators is demonstrated for fa = 0, 
02 = 7r/2 and 09 = — tt/2, and for e = 0.1. A graph on the 
right side of Fig. |2|is a graph of the Lyapunov function V (0). 
This function decreases but is not bounded in this example. 

B. Analytic Synchronization Condition for System 0 with 
Equal Coupling Strengths 

In this subsection we consider a special case of system 
([2]) when the coupling strengths are equal for all connected 
oscillators, i.e. we study the following system: 

4>i = h— y2 si n (4>j - 4>i)- (7) 

n jeN , 

The main result of this subsection is Theorem 1 which 
contains requirements on the initial phases and the coupling 
strength such that condition ([6| of Proposition 1 is satisfied 
and therefore system 0 achieves frequency synchronization. 



















Theorem 1 If D is a constant satisfying 0 < Dq < D < 
7 r, and 


K > 


n • \uji — ujj 

Pij • sin D 


( 8 ) 


for all i,j = l,...,n, then D t < D V t > 0 for the system 0 
in which the underlying topology is a graph with diameter < 2, 
and this system achieves frequency synchronization. 


Proof: Assume that at time moment T > 0, the value of 
D t is equal to D and before this moment it never exceeded 
D, i.e. D t < D \/t < T. We will show that under the 
conditions of this theorem, the maximum phase difference does 
not start to increase at time T by showing that requirement 
^ of Proposition 1 is satisfied. This will guarantee that the 
maximum phase difference D t will be always bounded by D. 

Condition must be satisfied for every two oscillators k 
and l with ^ = ^ nax and <f>J = 

tit ~ tit =^k~ ui 

~ ^ E sin (tik ~tit)~^E si n O j ~ tit) < 0- 

ieN k jeN t 


This condition will be satisfied if K > n T j UJk . and if we 

— PfcrsmD 

can show that: 


sin (^ “ + ^2 sin (^J “ > p ki • sinD. (9) 

ieN k jeNi 

Because and are respectively the maximum and mini¬ 
mum phase values at time T (see Fig. [3]): 

0 < - 4>J < D < 7T, 

o <fj -ff < D < IT, 

where 1 < i, j < n. Therefore, each summand in the left side 
of the inequality ([9]) is nonnegative. 

If vertices k and l are connected by an edge, both sums 
contain sin(</>^ — fj) = sinD, and thus the left side of ([9]) 
contains 2 sin IT 


Assume now that vertices k and l have a common neighbor 
- vertex m. Then, the left side of inequality contains the 
following sum: 

sin {tit ~ titn) + sin (titn ~ tit) 

= 2 S in(^).cos(^±4^) 

>2sin(^).cos(ti) 

= sin(^ - 4>f) = sin D. 

Inequality above holds because sin (^T)> 

0, and 

tit±tit_z ‘ttita 

2 2 2 — 2 2 5 

so that 

cos{ K±*Lz^ )> _ c< K^L >m{ ° ) , 


Therefore, the left side of ([9| contains a sum that is greater 
or equal than sin D for each common neighbor m of vertices 



Fig. 3: Oscillators and with the maximum and minimum 
phases, respectively 


k and /. In addition, if k and l are connected by the edge, there 
is a term 2 sin D in the left side of (|9|), and thus, inequality 

holds. This proves that condition ^6]) of Proposition 1 is 
satisfied under the theorem’s conditions. ■ 

Remark 1 If Do < |, the smallest value of boun d |8| ) will 
be achieved for D = When | < Dq <i r, bound takes 
its smallest value if D = D 0 . 

Remark 2 In the case of a complete graph, P Z] = n for 
each pair i, j of vertices, and the sufficient condition on K is 
the following: K > for all i, j. This bound coincides 

with the bound obtained in [8| for a complete graph. 

Remark 3 When diameter of a graph is larger than two, 
Theorem 1 cannot be applied in general, and condition ([6| of 
the Proposition 1 can be violated. For instance, if the distance 
between vertices k and j is more than two, then in condition ([6]) 
both sums may be equal to zero: sin — (j>\) = 0 Vi G 
sin((/>*• - 4>\) = 0 Vj G Ni, and - fj) > 0 if O k > 
uji. However, Theorem 1 can be applied to the graphs with a 
diameter more than two if every two oscillators with a shortest 
path between them of a length more than two, have equal 
frequencies. In this case condition ([6]) is always satisfied for 
such two oscillators. Indeed, if O k = Oi and (<^ — fj) = 

{timax -timin) = D < n, then ti k ~tii< 0 , because sin(0 t fe - 
4>\) > 0 and sin(<^- — fj) > 0 for all i G N k and j G TV/. 


C. Optimization Approach for System with Non-equal 
Coupling Strengths 

In this subsection the equal coupling strength assumption 
is relaxed. Instead of one coupling parameter K as was in the 
previous subsection, there are now \E\ coupling parameters 
Kij, where \E\ is the cardinality of the graph’s edge set E. 
Similarly to condition ^ in the Theorem 1, we will find 
bounds on the coupling strengths K tJ to guarantee frequency 
synchronization of system but instead of providing an 
analytic condition we will formulate an optimization 
problem whose solution contains the coupling strengths 
that guarantee <01 and are sufficient for synchronization. 

While in the Theorem 1 the goal was to find the minimum 
value of the coupling parameter K that guarantees synchro¬ 
nization, minimizing the sum of all coupling strengths ^ K tJ 

ijeE 













will be the goal for the case of non-equal coupling strength^] 

In condition ([6]) we assume that (\>\ = (frmax* = fimin 
and (f)^ — (j)\ = D. Since D < tt, all values of sin() functions 
in each sum of ([6]) are nonnegative. Instead of condition ([6]) 
we will consider a more strict condition on the coupling 
parameters, where we keep only summands corresponding to 
the neighbor oscillators of both oscillators k and Z: 

4>k-4>i <u k -Ui -^ • sin D 

~ E (•“■ • Sin (^fe “ #») + ^ • sin (0m - 0f)) < °> 

meN k i 

( 10 ) 

where e= 7V& D TV/ - is the set of common neighbors of 
oscillators k and Z. If there is no edge kl between oscillators 
k and Z, then K k i = 0 in ©• We will introduce constraints 
that do not contain phases and guarantee that condition (T0| ) 
(and ^ as well) is satisfied for all phase values. Optimization 
problem, whose \E\ variables are the coupling strengths Kij 
( ij G E ), that allowed to take nonnegative values, is formu¬ 
lated as follows: 


minimize: Kij, 

ijeE 

subject to: I Ok — OA — ^ kl . sinD 

n 

gm j J ^ / \ 

- - — [dm ’ Kkm + (1 — 8 m ) ’ Kim) E 0, 

meN k i 


(ID 


where 1 < fc, Z < n, and each S m may take values { 0 , 1 }. 
Since either S m or (1 — 8 m ) takes a zero value, variables K km 
and Kim do not appear together in each constraint. 


For each possible combination of values of S m there is a 
corresponding constraint, and, therefore for each pair of oscil¬ 
lators k and Z there are 2^*1 constraints in the optimization 
problem, where \N k i\ is the number of common neighbors of 
oscillators k and Z. For example, suppose that oscillators k 
and Z are connected and have a single common neighbor m, 
then optimization problem (TT| ) will contain two constraints 
for oscillators k and Z: 


I Wk -ui | - 
|^/c -Ul\ - 


2Kki 

n 

2 Kjd 

n 


• sin D — 


• sin D — 


sinD 

n 

sinD 

n 


Kkm < 0 , when 


Kim < 0 , when 


8m 

8m 



If, for example, oscillators k and Z are not connected and 
have two common neighbors mi, m 2 , then there will be four 
constraints for k and Z: 




sin D 


n 


{^K-km\ T~ Kkm2^j— 0? (8 t 


m\ — 8 m, 2 — 1/5 


\^k ^l\ ji ('^ m i Kim 2 ''j — O 5 (^mi — 8m ,2 — 0)- 

Thus, optimization problem GD contains in total 
2^*1 constraints. Although the number of constraints 

1 <k<l<n 

can be exponential in number of oscillators n, for some types 
of graphs it is polynomial in n. For example, for the graphs 
with star-tree topology, each pair of oscillators has at most one 
common neighbor, and thus, not more than two corresponding 
constraints. 

Remark 4 If all coupling strengths are required to be equal 
in optimization problem HU, then its solution is bound 
from the Theorem 1. Indeed, when all coupling strengths are 
equal, then Kki = Kkm = Kim in the constraint of GD for 

w k and Wi, J2 ($m ■ K km + (1 - S m ) ■ K tm \ = \N ki \, and 

meNki ' _ ' 

the constraint becomes: K > n _\ u:k ~ u }\ 

— P k i -sin D 

We will now show that solution to this optimization prob¬ 
lem satisfies conditions ( fTO] ) for all possible phase values. 

Theorem 2 Solution to the optimization problem GD 
satisfies conditions & and system achieves frequency 
synchronization. 

Proof: Suppose that K f , w here ij e E is a solution 
of the optimization problem fTT] ). We are going to show that 
condition fTO] ) is satisfied for two arbitrary oscillators k and Z 
with = D. This would imply that condition ( [ 6 t i s also 

satisfied since condition ( fTT)] ) is more restrictive than ^ 6 ]). 

For arbitrary phases <^, such that for all 

m G Nki, from condition 

- - 2 K ll • 71 

Wfe — u>i - — • sm D 

n 

~ E • sin (4 - fa) + ^ ■ sin(C - fa) 

meNki 

/ - - 2K kl • r, 

< uJk ~ ui -— • sm D 

n 

^ K?J (gin( ^ _ ^ + ^ 

meNki 

_ ^Kfoi . _ sin D • / t yr* \ 

<uj k ~ui - -smD -— ^ mm (K k m,K lrn ), 

n 11 meN kl 

because - </ 4 ) € [ 0 ,£>], {4> l m - $) € [ 0 ,£>], and - 

4L) + {4L ~ tf) = D. 

Now we can observe that for the right side of the last in¬ 
equality there exists a constraint in ( fTT] ) that guarantees that the 
right side is non-positive. If, for example, min(Ff^ m ,iT z * m ) = 
K^ m , then corresponding constraint in fTT] ) has 8 m = 1, 
otherwise 8 m = 0 . ■ 


\^k ^l\ ji K k m2 ^ — 0 ? (^mi — O 5 8 mi 2 — 1 ) •> 

sin D / \ 

\iC k CJi\ — yK k mi H - Kim 2 J — O 5 {8m 1 1 ? 8m 2 0 )j 


Objective function used in this article, therefore, is the F-norm of a vector 
of all coupling strengths. Other options could be employed, for example, 
the maximum norm Z°° could be used that corresponds to minimizing the 
maximum coupling strength. 


We finish this section with an example for which we 
found values of the coupling strengths that are sufficient for 
synchronization: first, under the condition that all coupling 
strengths must be equal and using the Theorem 1, and then, 
assuming that the strengths are allowed to be non-equal and 
solving the optimization problem GD- 

Example 2 In this example we consider four oscillators 
connected as shown on Fig. [4] and with following frequencies: 
























TABLE I: Synchronization conditions in our comparative analysis 



Bound on Coupling Strength 

Constraint on Initial Phases 

Analytic condition (Theorem 1) 

K — Pij -sin D (yi’j) 

Do < 7T 

Numerical condition (Theorem 2) 

Solution to (TT) 

D 0 < 7T 

Condition from | sj 

2 n- ||-B;T u>|| 

x 11 1,2 

A 2 ■7r-sinc(7 maa; ) 

Do < 7r; ||bJ</>(0)|| 2 < 7T 

Condition from |l4| 

X \ V / 2cr(u>) 

^ L* -sin D 

Do < tt; Z) 0° = 0; So < ^ 

i =1 

Numerical condition from [l2| 

K > 

yj £q ■ L- sin D 

X \ ™ , f 

— E sin (0fc-0i)+ E s in(cE 4>i) 

ie N k jeN l 

So <D 2 < 7T 2 


Oi — —0.1, O 2 = 0, CC 3 = 0.3, O 4 = —0.2. There are four 
edges in this graph, i.e. four coupling strengths Kij, and thus 
four variables in problem ( pTj ). Notice, that Pi 2 = Pi 3 = 
P 23 = 3, P 14 = 2 , and P 24 = P 34 = 1. If we assume that all 
the coupling strengths are equal, then by the Theorem 1 from 
previous subsection, sufficient for synchronization value of the 
coupling strength is: K = 0.5 • n = 2 (from the inequality for 
pair 34). Then, the sum of all coupling strengths is 4 • 2 = 8 . 

If we let the coupling strengths be different for the different 
edges, the optimization problem has a solution: K 12 = 0 . 8 , 
Kis = 2, K 2 3 = 0.2, and If 14 = 2. Now the sum of the 
coupling strengths is 5. Optimization problem for this example 
contains eleven inequality constraints (besides the constraints 
> 0 ). 

For optimization we used Matlab’s R2012a fmincon func¬ 
tion with default options. 

IV. Numerical Simulations 

In this section we present the results of simulations per¬ 
formed to demonstrate that for the graphs of diameter two, 
synchronization condition formulated in Theorem 1 is a less 
restrictive condition compared to existing ones. Since Theo¬ 
rem 1 guarantees existence of a Positively Invariant Set and 
frequency synchronization of system we compared our 
bound with the similar conditions that also guarantee existence 
of a PIS and frequency synchronization. To the best of our 
knowledge, there are three such conditions: Theorem 4.6 from 


1 4 

EV 

2 3 

Fig. 4: Connections between oscillators in Example 2 


(8), results from Q4), and conditions (analytic and numerical) 
in [12] . Therefore, we did not include into comparison analysis 
conditions from ® 0 and pi that only provide existence 
of an equilibrium and local stability. The numerical condition 
of p| is less restrictive then the analytic condition of the 
same article, and we here consider only the former one. In 
addition, we added to our comparison analysis a numerical 
synchronization condition from Theorem 2, which allows the 
coupling strengths to be different, and for each given example 
we calculated an average coupling strength of the solution to 



Each of the five synchronization conditions compared in 
this section consists of a bound on the coupling strength and 
constraints on the initial phases of oscillators. In particular, all 
synchronization conditions require that the difference between 
any two initial phases is less than 7 r (i.e. D 0 < 7 r). Addition¬ 
ally, synchronization conditions from |8), p2| and p4) have 
their own special constraints on the initial phases. The bounds 
on the coupling strength and corresponding requirements on 
the initial phases are summarized in Table 1. 

In the simulations we assigned a value of maxl^Do} 
to the constant D for our synchronization condition, because 
in this case bound © is the least restrictive as mentioned 
in Remark 1. In the bound from |8|, A 2 is the algebraic 
connectivity of a given graph, B c G M nxn ( n - 1 )/ 2 i s the 
incidence matrix of the complete graph with n nodes, O is 
a vector of frequencies, 0(0) - vector of initial phases and 
7 max = max{§, \\Bjm\\ 2 }- In the condition from |14|, £ 0 
is the squared Euclidean norm of a vector of the initial phases: 

i=1 

ct(uj) denoted he Euclidean norm of a vector of the intrinsic 
frequencies deviations: 
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(a) Fractions of random samples of initial phases that 
satisfy initial phase constraints 



(b) Average values of bounds in logarithmic scale 


Fig. 5: Simulation results of Experiments 1 and 2 


D is a constant whose value is defined as max{|, y/2£o}, and 
L* is defined as L* = 1 +diam ^Gy\ E c^ G ^ , where diam(G) = 2 
is the diameter of a graph G and \E C (G)\ is the cardinality of 
the set E C (G ) defined as 


E c (G)±E comp \E(G ), 

where E cornp is a set of n ^ n ~ 1 ^ edges of a complete graph 
with n nodes. 

In the condition from [12], cr(u)) and So are defined as in 
P4l , L = L *, and D = D 0 . 

In our analysis we compared the requirements on both, 
the initial phases, and on the coupling strength of the five 
synchronization conditions. 


Experiment 1 (comparison of the constraints on initial 
phases). In the first experiment we checked the restrictiveness 
of constraints on the initial phases of each of five synchroniza¬ 
tion conditions under consideration. We created 10 5 samples 
of the initial phases such that each phase was chosen from 
the (0,7r) interval. Then, for each sample we subtracted its 
mean phase value from each phase belonging to this sample. 
Therefore, the sum of the initial phases was equal to zero, and 
the maximum phase difference was less than n in each sample. 
We shifted the phase values of each sample by the sample’s 


mean because condition from [ 14 requires that ^2 $ = 0, 

and other synchronization conditions only depend on the 
relative values of the initial phases and thus are rotationally 
invariant. 


Next, for each sample we checked if it satisfies the con¬ 
straints on the initial phases of the synchronization conditions, 
and for each condition we calculated fractions of samples that 
satisfy its initial phase requirements. 

We repeated this experiment for different numbers of 
oscillators n in the system: n = 5,..., 10 and the experi¬ 
ment’s results are shown on Fig. [5a| Since our synchronization 
conditions in Theorems 1 and 2 do not contain any additional 


requirements on the initial phases, they can be applied for each 
generated sample of phases, and thus the fraction of acceptable 
initial phases is equal to one for all values of n. 

Fractions of acceptable initial phases for conditions from 
(14) and G3 monotonically decrease with the number of 
oscillators n as can be observed on Fig. |5a| 

Experiment 2 (comparison of the bounds on coupling 
strength). In the second experiment we compared the bounds 
on the coupling strength. For each fixed number of oscillators 
n = 5,..., 9 we randomly created 1000 graphs with n 
vertices and of diameter two. The initial edge set of each 
graph was empty, and we successively added random edges 
to it until the diameter was equal to two. For each graph 
we then created a random sample of initial phases, a random 
sample of frequencies, and calculated the bounds on K for 
each condition. For the numerical condition in the Theorem 2 
we calculated an average value of K for each example. The 
average values of bounds for each of five conditions under 
comparison are plotted on Fig. [5b] in logarithmic scale. In this 
experiment we sampled values of the frequencies from (0,1) 
interval, but the relative performance of the bounds does not 
noticeably change with the interval. 

The simulation results of Experiments 1 and 2 show that 
for graphs of diameter two our synchronization condition 
formulated in Theorem 1 is less restrictive in terms of both, 
initial phases and coupling strength compared to the existing 
conditions. Additionally, optimization-based condition in The¬ 
orem 2 provides a further improvement if the value of its bound 
is defined as the average coupling strength for each example. 

V. Conclusion 

In this article we employed the notion of a Positively 
Invariant Set to find a sufficient condition for frequency syn¬ 
chronization of heterogeneous Kuramoto oscillators connected 
by a graph of diameter two. We showed that an existence of a 
PIS ensures the boundedness of the trajectories of oscillators, 
which in turn, provides synchronization. For the case when 



























































the coupling strength is the same for every two connected 
oscillators, we provided an analytic synchronization condition, 
and demonstrated with simulations that this condition is signif¬ 
icantly less restrictive than existing ones. For the case when the 
coupling is allowed to take distinct values for different pairs 
of oscillators, we formulated an optimization problem whose 
solution - a set of coupling strengths - guarantees frequency 
synchronization. 
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